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Abstract 



A system of coupled two logistic maps, one periodic and the other chaotic, 
is studied. It is found that with the variation of the coupling strength, the sys- 
tem displays several curious features such as the appearance of quadrupling 
of period, occurrence of isolated period three attractor and the coexistence 
of the Hopf and pitchfork bifurcations. Possible applications and extensions 
are discussed. 
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The chaos in higher dimensional system is one of the focal subject of 
physics today. Along with the approach starting from modeling physical sys- 
tem with many degrees of freedom, there emerged a new approach, developed 
by Kaneko to couple many one dimensional maps [1-3] to study the behavior 
of the system as a whole. However, this model can only be applied to study 
the dynamics of a single medium such as the pattern formation in a fluid. 
What happens if two media border on each other? One may naturally leads 
to the model of coupled logistic maps with different strength parameters. 
Thus it is appropriate to inquire whether loosening the condition of strict 
identity might bring any new feature while keeping both maps to be logistic 
to hold the redundant controlling parameters minimal. Even two logistic 
maps coupled to each other may serve as the dynamical model of driven cou- 
pled oscillators. It has been found that two coupled identical maps possess 
several characteristic features which are typical of higher dimensional chaos 
[1,4]- 

In this letter, we report the results of numerical investigation on the sys- 
tem of two logistic maps with different strength parameters such that the 
one map lies in a period one stable attractor or a bifurcation point and the 
other in chaotic region when decoupled. Several new features previously 
unobserved are found. Most notable among them is the appearance of a 
period four cycle straight from the stable period one cycle. The other pe- 
culiar feature is the almost simultaneous occurrence of periods four, eight 
and sixteen right after the Hopf bifurcation which results in a very intriguing 
metamorphose of the attractor when one changes the coupling parameter. 

The system we study is two linearly coupled maps 

x n +i = (l-e)f(n,x n ) + ef(v,y n ) 

Vn+i = ef(n,x n ) + (l-e)f(y,y n ) (1) 

where the map / is taken to be the logistic map [5] with strength parameters 
fi. 

f(fi, x) = fix(l - x) (2) 

In case of /j, = v, two maps soon become synchronized no matter what ini- 
tial conditions may be, i.e., coupled maps are identical with a single logistic 
map. The interesting is the case c&n^v. In the following we fix the strength 
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parameters above and below critical value fi* = 3.56994 . . . for \i and v respec- 
tively. We choose the strength parameters fj, and v and regard the coupling 
parameter e as the controlling parameter. 

In Fig. 1, the attractors of the coupled- map are displayed as functions of 
coupling e. Fig. 1(a) shows the result of fi = 4 and v — 3, and (b) that of 
H — 4 and v — 2. They are two typical example of various values of \x and 
v. One immediately notices several interesting features. The fact that there 
are two chaotic regions in both e = and e = 1 ends seems odd at first sight, 
but after some reflection, one realizes that very weak e means very strong 
(1 — e) which brings chaos first to the variable x, then to y however weak the 
coupling term may be. The most salient feature is the appearance of stable 
period four cycle right after the period one around e = 0.77 in Fig. 1(a). 
Another, found in both (a) and (b) cases, is the sudden filling of the x and y 
space around e = 0.85 and above. The broad window-like region with period 
four around e = 0.9 in the case of (b) is also noteworthy. 

Identifying stable and unstable periodic points often gives the skeleton 
view of the dynamics of the system. With the vector notation of the variables 
z— (x, y) and the operator notation of the map 

z n+ i = F(z n ), (3) 
one can define the periodic points of cycle N as the solution of the equation 

z N = F N (z N ), (4) 

N times 

Z * 

where F N (z N ) is defined by F N (z N ) = F(F(F ■■■F(z N )-- •))• Bifurcations 
occur at the parameter value of e where the periodic point satisfies 



E, 



max 



dF N _ 



1, (5) 



in addition to eq.(4). Here, E max [G] is the maximum eigenvalue of the ma- 
trix G. We numerically solve eq.(4) for iV = 3,4, and 8 for entire range of 
e = ~ 1. The results are presented in Fig. 2 where the unstable orbits are 
shown as well as the stable ones. One can clearly see from iV = 4 and N = 8 
skeleton, that the period four occurring at e = 0.77 is the result of the folded 
period doubling bifurcation sequence whose period two starts at e = 0.75 after 
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the quadrupling of the period when one moves from period one at larger e 
to a higher periodic orbits at smaller e. While this type of inversion of the 
sequence is nothing of unimaginable nature in principle, it has never been 
observed in simple maps to our knowledge. One recognizes easily that it is 
rather hard to obtain in a usual one-dimensional map since it requires spe- 
cial tuning of functional dependence of the map on the controlling parameter 
to ensure that the condition eq.(4) is satisfied for N = 4 before N = 2. We 
therefore think it significant that it emerges from nothing but linear coupling 
of two logistic maps of different strength parameters. One might presume 
that it should be observed rather frequently in higher dimensional coupled 
lattice map models once the requirement of strict uniformity of the elements 
is lifted. In fact, preliminary investigation on three coupled logistic maps 
supports this view. 

There are two notable facts in the periodic cycle diagram above e = 0.85. 
One is that no periodic point exists in the region e < 0.86. The examination 
of the x — y profile shows that the onset of "filling area" at approximately 
e = 0.85 is the result of the Hopf bifurcation [6]. The second is that around 
e = 0.88, period four and period eight start (as well as other higher periods 
not shown here) appearing almost simultaneously. While these occurrences 
fit into the generic transition scenario of "cycle one — > Hopf — > cycles — >■ . . . 
— >• chaos" found by Kaneko in an early study [1,3], the second fact, crowded 
onset of many different periodic cycles makes this transition to chaos a very 
intricate one. These points are visually displayed in the phase portrait of the 
attractor at several values at and above e = 0.85 shown in Fig. 3. The first 
two figures (a) and (b) clearly show the Hopf bifurcation around e = 0.85. 
Subsequent distortion and transformation to an aurora-like strange attractor 
observed in Fig. 3(c) to (f) certainly match those found in far more com- 
plicated system in its aesthetic appeal. One curious feature found both in 
Figs. 2 and 3(b) is the early appearance of period 3 cycle around e = 0.87. 
At present, we are unable to notice any direct role of this period to the shape 
and stability of the attractors around this region [7]. 

In summary, we have constructed a system consisting of chaotic and pe- 
riodic maps coupled together. While the dynamics of the system does not 
go beyond the known bifurcation schemes in the two-dimensional dissipative 
system (which is certainly not to be expected), it is found that various bi- 
furcations occur in such a combination to give the system several intriguing 
features. 
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Finally, few words on possible applications and extensions are in order. 
It should not be too difficult to construct a circuit (electronic, for example) 
to materialize the system proposed here. Actual applicational usage of qua- 
drupling of stable states might be envisioned. Also, replacing the logistic 
map with other maps, cycle map or quadratic map for example, might be 
interesting to see the generality and/or the new aspects of the findings here. 
Another application is possible by increasing the number of maps coupled 
each other. The coupled lattice map having "impurity" , or one element with 
different strength parameter from all the rest might exhibit non-conventional 
features. A lattice with alternating strength might also be an interesting sys- 
tem. In a word, loosening the condition of strict uniformity of the elements 
might bing some new feature to the coupled lattice map. 

We express our gratitude to O. Morimatsu for the supports in computa- 
tional aspects of the work. The numerical calculation has been performed on 
VAX3100 of the Theory Division, Institute for Nuclear Study, VAX6000/400 
at Nuclear Science Center, Department of Physics, University of Tokyo and 
also on Sparc Station 2 at Department of Physics, Hosei University. 
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Figure Captions 



Fig. 1 The phase diagram of the coupled maps eqs.(l)-(2). (a) is for the 
parameter values n = 4.0, v = 3.0, and (b) is for mu = 4.0, v = 2.0. 

Fig. 2 The "skelton" of Fig. 1, showing the position of the periodic cycles of 
the maps, (a) is for /i = 4.0, v = 3.0, and (b) for /i = 4.0, v = 2.0 as 
before. The caption number indicates the period of the cycle. 

Fig. 3 Spatial portraits of the attractors for the coupled logistic map with 
H = 4.0 and v = 3.0 at e = 0.852 (a), e = 0.868 (b), e = 0.8757 (c), 
e = 0.877 (d), e = 0.885 (e), and e = 0.900 (f). 
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